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ABSTRACT 

A numerical simulation model of the injection and cooling of nonthermal particles energized by 
shocks formed in merging clusters of galaxies is used to fit radio and X-ray data observed from the 
Coma cluster of galaxies. The results are consistent with a primary merger-shock origin for both 
the diffuse radio halo emission and the hard X-ray excess measured with Beppo-SAX and RXTE. 
For equal (1%) efficiency of power injected in nonthermal protons and electrons, we predict that the 
Coma cluster of galaxies will be significantly detected with the space-based observatory GLAST, and 
marginally detectable with the ground-based 7-ray observatories VERITAS and HESS. Significant 
TeV detections are possible if the nonthermal proton intensity is greater due to a larger efficiency of 
nonthermal hadron acceleration, or to past merger events. The nonthermal hadronic content in Coma 
is also revealed by a weak, hard secondary emission component at ~ 10 - 100 GHz. The spectral 
softening of the radio emission at large radii from the Coma cluster core derives in this scenario 
from the decreasing shock speed away from cluster center for an on-axis merger event. We discuss 
differences between merger and delayed turbulence models. 
Subject headings: 



1. INTRODUCTION 

The discovery of the first radio halo, Coma C, demon- 
strated the existence of a population of nonthermal elec- 
trons within the intra cluster medium (IC M) of the Coma 
cluster of galaxies fsee lSarazinl l)1988f > and lBiviano . ( 1998) 
for review; rad io intensity maps and spectra o f Com a 
are reported bvlThierbach. Klein, fc Wielebinskil iI2003|L 
ISchlickeiser et alJ lj!987|) . and IGiovannini et alJ (|1993|) i 
The origin of these energet ic electrons co uld be associ- 
ated with a merger event l)Tribble|ll993|) . for example, 
between the central NGC 4889 and NGC 4874 galaxy 
groups. In addition, the radio relic 1253+275 is found 
on the Coma cluster (A1656) periphery, possibly indi- 
cating the presence of a shock resulting from a recent 
or ongoing merger event. Several tens of galaxy clus- 
ters, representing ~ 10% of the cluster population, are 
reported to have radio halos, and th ese clusters p refer- 
cntially show recent merger activity l)Buotd l2001). Be- 
cause of their relatively short radiation lifetimes, GHz- 
emitting nontherm al electrons have to be accelera ted on 
timescales shorter ijSarazi n1 li999t |Pet7osian1l200l than 
merger timescales. 

Optical measurements of radial velocities of galaxies 
reveal the cluster dynamics dominated by the gravita- 
tional potential well of the dark-matter halo. Although 
Coma is considered to be the nearest rich relaxed cluster, 
with a FWHM diffuse X-ray angular extent ~ 30', opti- 
cal surveys of the radial velocities of 1174 galaxies in the 
Coma cluster llColless fc Dunnlll99fitlBiviano et al]ll99ft 
lEdwards et al.ll2002fl show that, it is in a highly dynamic 
state as a result of a recent merg er event. 

Observations with ROSAT ^ White. Briel. fe Henrvl 
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1993) and ASCA ijHonda et alJ (|1996j) : see references in 
ijBriel et al.N200lfO show multiple X-ray emission peaks 
from the core of the Coma cluster. The diffuse clus- 
ter emission peak associated with the cluster center is 
coincident with the galaxy NGC 4874. A second diffuse 
emission peak was found approximately halfway between 
NGC 4874 and NGC 4889. XMM-N ewton observations 
ijArnaud et alJ200HlBriel et al.l200l|) did not detect this 
second peak but found X-ray emission peaks coincident 
with both NGC 4874 and NGC 4889, a relaxed ther- 
mal structure within 10' of the cluster center, multiple 
point sources, and stripped gas on the cluster periphery. 
Chandra observations of galaxy clusters with radio halos 
reveal complex temper ature maps indicat ing recent or 
ongoing merger events ijGovoni et alJl2004j) . In the case 
of Coma, Chandra results indicate that the X-ray emis- 
sion from the coronae surrounding the dominant galaxies 
NGC 4874 and NGC 4889 and their galaxy groups arises 
from denser, X-ray gas clouds with temperatures of a 
few keV that are in eq uilibrium with the intracluster gas 
ijVikhlinin et alJl200l. 

BevvoSAX l|Fusco-Femiano et all |2004|) and Rossi 
X-Ra y Timing Explorer {RXTE) ijRephaeli fe Gruberl 
2002) observations of the Coma Cluster show hard X-ray 
(HXR) emission in excess of single-temperature brems- 
strahlu ng, though its detection with BeppoSAX is dis- 
puted (Rossetti & Molcndi 2004), and the hard tail re- 
mains unconfirmed with INTEGRAL. BeppoSAX results 
ijNevalainen et alJ2004|) also indicate that clusters under- 
going recent or ongoing merger events have HXR excesses 
with 20 - 80 keV luminosities « 10 43 - 10 44 ergs, whereas 
fully relaxed clusters are consistent with no HXR excess. 

At 7-ray energies , the Coma cluster was not detected 
ijReimer et al.ll2003|) with the EGRET telescope on the 
Compton Observatory at > 100 MeV. Statistical analy- 
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sis of the positional association of unresolved gamma- 
ray emission from rich clusters of galaxies in the lo- 
cal universe shows a weak (~ 3ct) positive correlation 
( Sch arf fe Mukheriee^OO^ . In addition, there are claims 
of an association between unidentified, high galactic lat- 
itude (|6| > 20°) EGRET so u rces and galaxy clus - 
ters l|Kawasaki fc Totanil J200l: iColafr ancesco (2002), 
though contested bv lReimer et alJ lf2003j) ). The results of 
searches for TeV emission with the ground-based imaging 
air Cherenkov telescopes (IACTs) HESS (High Energy 
Stereoscopic System) and VERITAS (Very Energetic Ra- 
diation Imaging Telescope Array System) as they reach 
their design capabilities are anxiously awaited, 1 but the 
sensitivities of IACTs and coded mask telescopes such as 
INTEGRAL are compromised by the extent of Coma's 
diffuse emission iJGabici & Blasi 2004) which, as in the 
case of the thermal X-rays, could subtend a large frac- 
tion of a square degree. The diffuse extent of the non- 
thermal emission in clusters of galaxies, as indicated by 
radio halos and relics, also makes it more difficult to map 
emission from cluster and accretion shocks with narrow 
field-of-view pointing telescopes. The > TeV 7-rays from 
Coma will also be attenuated by the diffu se extragalactic 
infrared radiation ijGabici fc Blasill2004[) . 

There have been numerous attempts to model the dif- 
fuse radio emission associated with clust er radio halos 
in the context of a cluster merger model. iMiniati et alJ 
l)2001|) treat acceleration of cosmic-ray protons in the 
context of a cosmological structure-formation simula- 
tion. They find that strong shocks form from cool, 
low density ICM gas accreting onto the cluster periph- 
ery. Such shocks are not expected to penetrate into 
the central regi on of the cluster environment. Numer- 
ical s imulations llRoettiger. Burns, fc Stondll99?l iRickerJ 
1998; Rickcr fc Sarazinll2001j) show that shocks forming 
in merger events between two virialized clusters of galax- 
ie s can, however, penet rate into the cluster core. 

iGabici fc Blasil l)2003) found that shocks resulting from 
a major merger event (as defined in Section 2) between 
two virialized clusters of similar mass produce nonther- 
mal particle spectra that are too soft to account for the 
observed emission from most radio halos. Minor merg- 
ers are more effective for injecting hard electron s pectra. 
Berrington & Dermer (2003) (hereafter |BD| l|2003in show 
that the typical peak Mach speeds M. of cluster merger 
shocks in minor mergers are M ss 3 - 5. Primary elec- 
tron synchrotron radiation injected by merger shocks can 
account for the observed radio emission from Coma C, 
as we show here. 

Besides a directly accelerated electron population from 
a merger shock, nonthermal radio halos detected from 
Coma and more than ~ 20 clusters of galaxies could also 
originate from the decay produc ts of proton interactions 
|DennisorJll98a IVestrandlll987|) . IBlasi fc Colafrancescol 
( 1999) performed a detailed study of this model for the 
Coma cluster and find that if the radio emission is due to 
secondary electrons, then EGRET should have detected 
Coma. Moreover, this model is difficult to reconcile with 
the correlation o f cluster radio halos with merger activ- 
ity l)Buotell2001[) or the softening of Coma's radio spec- 

1 The latitudes of Whipple/ VERITAS and HESS are +31.4° and 
—23.3°, respectively, so Coma, at declination +38.0°, is a better 
target for VERITAS than HESS. 



trum with radius, because clusters of galaxies will con- 
fine protons for the age of the universe l|Volk et allll996t 
iBerezinskv et alJ ll"9971. The cluster radio halos could 
also be a consequence of delayed acceleration by second- 
order Fermi proc esses associated with m agnetic turbu- 
lence in the ICM iSchlickeiser et alJl9 87). The most de- 
tailed study to date l)Brunetti*et alJ l2004') of the plasma 
physical processes results in a model for the Coma radio 
spectrum that could fit the data. Some concerns with 
the merger model, and fi ts to Coma data in a simplified 
merger model is given by Blasi (2001) (with references to 
earlier work). We consider these points in our Discussion 
section. 

In Section 2, we adapt our merger-shock model (iBDl 
2003) to the Coma cluster environment. Our approach 
directly ties the substructure in Coma to the radio and 
HXR excess. Assuming a merger-shock origin for the 
nonthermal radiation from the Coma cluster, we show in 
Section 3 that if protons are accelerated with the same 
efficiency as electrons, Coma is predicted to be detected 
at GeV energies with the Gamma Ray Large Area Space 
Telescope ( GLAST), and at TeV energies with the IACTs 
VERITAS and HESS. Hadrons rema ining from previous 
merger shocks l|Gabici fc Bl asi 2003) could increase the 
7-ray flux, but loss of sensitivity to extended sources 
could decrease detectability. An estimate of GLAST sen- 
sitivity for the Coma cluster in the Appendix shows that 
even for the most optimistic fluxes of 7-rays allowed by 
radio observations, Coma will be detected with only ~ 5a 
significance. 

The results are summarized in Section 4, where we 
discuss and compare predictions of merger and reaccel- 
eration models. We argue that gamma-ray detection of 
Coma will support a merger-shock origin of the nonther- 
mal radiation from Coma, though second-order turbu- 
lent re-acceleration effects could play an associated role 
in nonthermal emissions from clusters of galaxies. The 
Appendix gives some analytic results for cluster merger 
physics that provide a check on the numerical results. 

2. THE COMA CLUSTER ENVIRONMENT 

Because nonthermal spectral power depends sensitively 
on the parameters of the c luster envi ronment, we modi- 
fied the code developed bv lBDl l)2003fl specifically to the 
properties of Coma (see Table 1). The Coma cluster of 
galaxies is well-described by an event w here a dominant 
cluste r of total mass Mi = 0.8 x 10 15 M Q (jColless fc Dunnl 
1996) merges with a s maller cluster of total mass M? sa 
0.1 x 10 15 M^ c luster (jVikhlinin. Forman. fc Joneslll997t 
IDonnellv et afll99"l . We assume a gas mass fraction of 
«5%. 

Let A m represent the mass ratio between the merg- 
ing and dominant clusters. Merger events with mass 
ratio A m < 0.1 are designated acc retion events 
()Salvador-Sole. Solanes. fc Manriaueill998l) . The core of 
the domina nt cluster is found to r emain intact unless 
A m > 0.6 (jFuiita fc Sarazinll200H) . In our proposed 
model, the mass ratio of the two clusters is A m = 0.13. 
This places our model well below the threshold when the 
internal structure of the dominant core is destroyed, and 
classifies our model as a cluster merger event. In addi- 
tion, the mass ratio is great enough that the shock result- 
ing from th e merger event will penetrate into the core of 
the cluster (jRoettiger. Burns, fc Stonelll999j) . Our semi- 
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TABLE 1 
Properties of Merging Clusters 
in Coma 



Redshift z 


0.023 


d L (cm) 


3.0 x 10 26 cm 


Mi 


0.8 x 10 15 iW Q 


M 2 


0.1 x 10 15 M 


(Txh 


8.21 keV 


(Txh 


2.0 keV 


r c ,i 


0.257 Mpc 


r c ,2 


0.15 Mpc 


^*max,l 


1.3 Mpc 




0.75 Mpc 


P 


0.705 


P0,1 


6.39 x 10 -27 gm cm -3 


P0,2 


1.67 x 10~ 27 gm cm" 3 



analytical dynamics approach compares favorably with 
N-body simulation results until the centers-o f-mass of 
the clusters pass through each other (BD 2003). At later 
times, the dynamics of the merging cold dark matter ha- 
los and the hydrodynamics of the gas are not accurately 
treated, and during this phase we expect rapid coales- 
cence of the cold dark matter halo and shock quenching 
as the gas dynamics enters a Sedov phase. 

The calculation of nonthermal electron injection, elec- 
tron energy losses and the production of the Compton- 
scattered emission for the various radiation fields, in- 
cluding the cosmic microwave background radiation 
(CMBR), the stellar radiation, a nd the X-r ay radia- 
tion fiel d s, foll o ws the treatment of[BD (2003J (see also 
ISarazinl £l999); P etrosiar] lj2001ft ). For a galaxy lumi- 
nosity function described by the Schechter luminosity 
function with parameters M* = —21.26, = 107, and 
a = 1.25 l|Schechterll976ft . the energy density of the stel- 
lar radiation field is found to be 



4.87 x 10" 



R 



1 Mpc / 



MeV cm 



(1) 



where R is the radius of the Coma cluster. 

The Coma cluster has a tenuous ICM that emits ther- 
mal X-ray bremsstrahlung. To model the emission from 
the ICM, we assume that the ICM is well approximated 
by a (3 model 2 with core radius r c = 0.257 Mpc, central 
electron density p e o — 3.82 x 10~ 3 cm' 3 , central proton 
density po.i = 6.39 x 10 -27 gm cm -3 , power-law slope 
= 0.705, and a mean gas tempera ture (Tx) = 8.21 keV 
ijMohr. Mathiesen. &; Evrardl 11999ft (for the total mass 
modeling, we assume for simplicity that the gas is en- 
tirely composed of H). The energy density of the X-ray 
radiation field is then given by 



Px 



: 1.30 x 10~ 10 



— V 

1 Mpc J 



MeV. 



(2) 



This energy density corresponds to a total free-free lu- 
minosity of ~ 1 x 10 45 ergs s" 1 emitted within a cluster 
radius of 1 Mpc. The energy density of the CMBR is 



Pcmb = 2.7 x 10" 7 MeV cm" 



(3) 



2 As shown in Appendix A. 3, the use of an NFW profile gives 
relative merger speeds that are similar to that obtained with a 
/3-model. 



Given the evidence for a recent merger event in Coma, 
as described in Section ^ we argue that merger shocks 
accelerate nonthermal particles to produce nonthermal 
radio and X-ray emission. The evolution and proper- 
ties of the forward and reverse cluster merger shocks are 
determined by the relative speeds and densities of the 
merging clusters. We have calculated the relative clus- 
ter speeds and shock properties following the method de- 
scribed in BD (2003). The particle number density of the 
ICM is assumed to be well approximated by a spherically- 
symmetric model whose parameters are given above. 
The dark matter halo, whose total mass is normalized to 
Mi out to a maximum radius r max as given in Table 1, 
is assumed to follow a pr ofile similar to the cluster ICM. 
Beyy o-SAX observations ijEttori. De Grandi. fc Molendil 
12002ft support the claim that t he da rk matter halo of 
Coma is better fit by the lKind l)1962ft approximation to 
the isothermal model, given by 



Po 



-3/2 



(4) 



where po — CPo.i is the central density, £ is the ra- 
tio of the total matter density to the normal mat- 
ter density pn i and r r is the core radius, than the 
iNavarro. Frenk. & Whitel l)1997ft (NFW) density profile. 

The efficiency parameters rj e and r\ v giving the fraction 
of power dissipated by the shock in the form of nonther- 
mal electrons and protons, respectively, are set equal to 
1%. The injected particle spectrum is described by a 
power law in momentum space truncated by an expo- 
nential cutoff. The cutoff energy is determined by time 
available for acceleration since the onset of the shock, 
comparison of the size scale of the system with the par- 
ticle Larmor radius, and competition of the accelera- 
tion time scale with radiative-loss time scales (§2.2 in 
IBDI ((2003)). The total energy liberated from the cluster 
merger event and deposited into the nonthermal particles 
is ~ 1.5 x 10 61 ergs. 

3. MODEL RESULTS 

The model that best represents the multiwavelength 
data of the Coma Cluster is found to be described by the 
merger of two clusters nearing collision time, i C oiij when 
the centers of mass of the two cluster s are almost co- 
incid e nt, as supported by observations iColless h Dunnl 
fl996t iBiviano et all 11999: lEdwards et all 12002ft . With 
the masses of the merging clusters given above, and an 
initial separation of 2.25 Mpc, the observations take place 
at t co u w 1.0 x 10 9 yrs after the onset of the merger. The 
observed redshift of the Coma cluster is z a = 0.023. For 
this value of i co ii, the redshift of the cluster at the start 
of the cluster merger event (when the outer radii of the 
two /3-model profiles first intersect) for a flat ACDM cos- 
mology with fig = 0.3 and J7a = 0.7 is Zi — 0.10. Here 
we use a Hubble constant Hq = 75 km s' 1 Mpc' 1 , giving 
a luminosity distance = 3.0 x 10 26 cm for Coma. 

The total energy spectra of primary protons and pri- 
mary and secondary electrons (including positrons) that 
fit measurements of the Coma spectrum are shown in 
Fig. ^ As can be seen, the proton distribution is well 
described by a power-law spectrum N(E) cx E~ 2 2 be- 
tween 1 GeV and » 10 17 eV. At the highest energies, 
the spectrum is cut off due to the available time required 
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Fig. 1. — Total energy spectra of primary protons and primary 
and secondary electrons, plotted in the form of E 2 N(E), used to 
fit the measured spectral energy distribution of the Coma cluster. 
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to accelerate protons to the highest energies. The strong 
radiative losses, principally due to Compton scattering 
with the CMBR, causes the abrupt cutoff of the primary- 
electron spectrum. Because of the production kinematics 
and pile-up due to radiative losses, the secondary elec- 
tron spectrum displays a bump at a few GeV. 

Fig. |3 shows the different spectral components used to 
fit the multiwavelength Coma spectrum. For the single 
merger event calculated here, the primary electron radi- 
ation generally dominates secondary emissions at radio 
and X-ray wavelengths, but secondary emissions make 
a contribution at mm radio wavelengths and at > 100 
MeV energies for the assumed 1% efficiency fractions r] e 
and i] p of energy injected in the form of nonthermal par- 
ticles. The value of rj e is determined by normalizing the 
primary electron Compton emissions to the > 40 keV 
nonthermal X-ray flux, insofar as only a small fraction of 
the X-ray flux is produced by bremsstrahlung and Comp- 
ton emissions from secondary — ► processes. The 
hardening at « 1 GeV in the primary Compton compo- 
nent is stellar radiation photons scattered by the primary 
electrons. The ratio of the energy densities, and there- 
fore peak fluxes, is « 50, from eqs. and ©. Scatter- 
ing of the thermal diffuse X-rays to > GeV energies is 
suppressed by the absence of primary electrons at these 
energies (Fig. 0). 

At > 100 MeV 7-ray energies, most of the emission is 
from secondary 7r°-decay 7 rays formed by nuclear pro- 
cesses in collisions of nonthermal protons with the diffuse 
gas. The flux of the secondary emissions is in proportion 
to the assumed fraction, 1%, of the power swept-up by 
the forward and reverse shocks and injected in the form 
of a nonthermal power-law momentum distributions of 
hadrons. The total spectral energy distribution used to 
fit the Coma spectrum is shown in Fig. |3| where we vary 
the fraction r\ p of energy deposited in the form of non- 
thermal hadrons from 1% to r/ p = 10%. 

A comparison of the merger model wi th the measured 
ijThierbach. Klein, fc Wielebinskil 12003^ radio emission 
from the radio halo of Coma is shown in Figure 0] The 
mean magnetic field, obtained by normalizing to the ra- 
dio emission, is 0.22 /xG, and is in agreement with ana- 
lytic estimates (see Appendix). Our model is in accord 



Fig. 2. — Spectral emission components for the fit to the volume- 
integrated multiwavelength spectral energy distribution observed 
from the Coma cluster of galaxies. Synchrotron, Compton, and 
bremsstrahlung emission components are shown for the directly- 
accelerated primary nonthermal electrons (solid curves) and for 
nonthermal leptons and 7r°-decay 7 rays formed as secondaries in 
inelastic nuclear pion-producing processes (dotted curves). Also 
shown is the diffuse thermal bremsstrahlung emission (dashed 
curve) . 
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Fig. 3. — Total spectral energy distribution used to fit obser- 
vations of the Coma cluster with 1% of the power swept up by 
the forward and reverse shocks dissipated in nonthermal electrons. 
The two curves with r\ v = 1% and r\ v = 10% correspond to the 
fraction of shock power deposited in nonthermal protons. Separate 
components for the total primary (solid curve) and total secondary 
(dotted curve) emissions for the case r) p = 1% are also shown. 



with a primary electron source for the radio emission — 
a comparable contribution to the radio flux from sec- 
ondary nuclear production would require that rj p ap- 
proach 100% to produce the measured radio emission, 
and would furthermore require a very different nonther- 
mal proton spectrum. From spectral considerations, ob- 
servations of Coma imply that rj p < 20%; otherwise a 
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Frequency [MHz] 

Fig. 4. — Comparison of radio observations of the Coma 
cluster with results of the cluster merger-shock model. The 
solid curve is the total radio emission from the model, the 
dashed curve is the contribution to the total emission from pri- 
mary electrons, and the dotted curve from secondary electrons 
and positrons. The solid circles are observational data points 
ll'hierbach. Klein, fc VVielebinskil2003) . 



hard hadronic component would be seen in the multi- 
GHz radio emission. We note that the fluxes from sec- 
ondary electrons and positrons calculated here may over- 
estimate the actual fluxes from secondary leptons by as 
much as a factor of « 2 because a uniform density pro- 
file, appropriate to the central region of the cluster, is 
used in the calculation of secondary production. This is 
nevertheless a good approximation because the density 
in the central region is roughly constant and most of the 
nonthermal radiation is formed by interactions in this 
region. 

Our model assumes a single merger event. Because 
galaxy clusters are b uilt up by multiple merger events 
ijGabici fc" Blasi 2003), and GeV - PeV protons and ions 
do not escap e fro m galaxy clusters during a Hubble time 
fVolk et alJIlQQfit iBerezinskv et al] ll997). a nonthermal 
population of protons accelerated in previous merger 
events is likely. Coma may have seen ~ 10 merger events 
during its buildup. Previous merger events are less ener- 
getic and nonthermal protons experience collisional and 
catastrophic losses through interactions with the ther- 
mal protons in the ICM in the meantime, so the total 
energy enhancement of hadrons over leptons is probably 
not more that f=a 2 — 10 from past merger events. The 
superposition of proton spectra from th ese past events 
may also produce non-power law spectra ijGabici fc Blasil 
2003). This would also hold for the varying shock com- 
pression ratios in a single merger event, though we find 
that because most of the nonthermal hadrons are injected 
when the spectrum is hardest, the proton spectrum is 
accurately described by a power-law (Fig. Q with index 
close to the hardest reached during the merger event. As 
demonstrated in the Appendix, injection indices reaching 
ps 2.2 are feasible during a minor merger, 

The calculated thermal and non-thermal X-ray fluxes 
are compared in Fig. [5] with the HXR flux mea- 




Photon Energy (keV) 



Fig. 5. — Comparison of X-ray observations of the Coma cluster 
with the cluster merger-shock model. Data points from different 
instruments are labeled in the diagram. The OSSE data points are 
2<r upper limits. The solid curve is the sum of nonthermal emission 
(dotted curve), dominated by Compton-scattered CMB radiation 
from primary electrons, and thermal ((Tx) = 8.21 keV) brems- 
strahlung emission (short-dashed curve). Point-source sensitivities 
for a 3<r detection in 10 6 seconds with INTEGRAL are plotted, 
though the actual sensitivity is a factor 2 — 3 worse because of 
Coma's extent. 

sured from the central (~ 2.2 Mpc) region of the 
Coma cluster. Data points show the HXR flux mea- 
sured with the Phoswich Detect ion System on Beppo- 
SAX l|Fusco-Femiano et al.ll200l . the HEXTE (High- 
Energy X-Ray Timing Exper iment) instrument on 
RXTE iRephaeli fc Gruberll200l. and the OS SE 2a up- 
per limits (RcDhacli. Ulmer. fc Gruberl 11994). Despite 
the difference in spectral range and field-of-view of the 
Proportional Counter Array on Beppo-SAX and HEXTE 
on RXTE, both instruments indicate the presence of an 
HXR excess observed in Coma. The existence of the 
HXR ta il in the Coma cluster has , however, been chal- 
lenged (Ro ssetti fc Molendl l2004j) . Observations with 
ISGRI on INTEGRAL could help to resolve this con- 
troversy, as indicated by the sensitivity curves plotted 
in Fig. [3 But because the Coma Cluster is extended, 
the sensitivity is degraded by a factor of 2 - 3 by the 
use of a coded-mask imager on the INTEGRAL satellite 
(A. Vikhlinin, private communiation, 2004). Astro E-2 
(Suzaku) should resolve this important question in the 
near future. 

Fig. HO shows the predicted 7-ray emission from the 
Com a cluster of galaxie s. Sensitivity limits are taken 
from lWeekes et al.l pOO^ for the VERITAS point-source 
sensitivty (which is comparable to the HESS sensitiv- 
ity). Because of Coma's extent, the sensi tivities of these 
IACT s can be ps 4x worse at 1 TeV l|Gabici fc Blasil 
120041) than shown. The predicted 7-ray emission falls 
below the EGRET sensitivity curve and the measured 
2cr upper limit of 3. 81 x 10" 8 ph(> 100 MeV) cm" 2 s" 1 
(|Reimer et al.ll2003|) . Our model predicts that GLAST 
will significantly detect the non-thermal 7-rays from 
Coma to energies of several GeV. Reduction in the sensi- 
tivity of GLAST due to the Coma's extended structure is 
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Fig. 6. — Predicted 7-ray emission from the Coma cluster of 
galaxies from the cluster merger shock model. The solid curve is 
the predicted photon flux (in units of ph(> i? pca k) cm -2 s _1 ). 
Sensitivity limits for EGRET, MAGIC, GLAST, and VERITAS 
and HESS are included. The EGRET limits are for 2 weeks in 
the pointing mode, and the GLAST limits are for 1 year in its 
scanning mode. The quoted VERITAS, MAGIC an d HESS point- 
sourcc sensitivities are for 50 hour, 5<j observations IWeekes et alj 
2002), and will be degraded due to the angular extent of the Coma 
cluster's emission. 

not significant, as shown in the Appendix. Furthermore, 
we predict that both VERITAS and HESS could have 
high confidence (> 5a) detections if r] v = 10% or there is 
significant nonthermal protons left over from past merger 
events. Attenuation of sub-TeV 7-rays by the extragalac- 
tic diffuse infrared and optical radiation fields is not sig- 
nifica nt, though it is for multi-TeV fluxes ijGabici fc Blasil 
l200l . 

4. DISCUSSION 

We have developed a cluster merger model based on 
observations of merger activity in the central regions 
of the Coma cluster of galaxies. Our model employs a 
semi-analytic treatm ent of clus ter dynamics based on an 
N-body simulation ijBDI I2003|) where nonthermal elec- 
trons and protons are injected with an intensity and 
spectrum that depends on the strength of the evolving 
shocks that are formed by the cluster merger. This rep- 
resents a significant i mprovement over previous analytic 
merger models (e.g.. iBlasi &: Colafrancescol 119991 iBlasil 
I200H IReimer et al.l I2004|) bv considering the evolution 
of the particle distribution and system parameters with 
time. The model also provides a joint treatment of both 
electron and hadron populations and detailed spectral 
fits in comparison with those yet available with numeri- 
cal hy drod ynamic codes, such as given bv lMiniati et all 
( 2001) and lMiniatl l)2003|) . Moreover, our model is specif- 
ically tied to the observed substructure of Coma. We find 
that a detailed merger model gives a good fit to the mul- 
tifrequency spectrum of the Coma cluster, confirming the 
simpler analytic models. 

The measured spectrum of the diffuse radio emission 
of Coma C is consistent with a model of synchrotron ra- 
diation from electrons accelerated at merger shocks. Be- 



cause of the short radiative lifetimes (« 50 - 200 Myrs, 
depending on magnetic field) for GHz-emitting electrons, 
compared to the ^0.5 Gyr merger time scale (see Ap- 
pendix), the primary synchrotron emission emitted over 
the lifetime of the injected electrons is very soft, thus 
accounting for the soft radio spectrum shown in Fig. 0] 

The radio emission from Coma C shows an increasingly 
soft spectral index wit h radius ijSchhckeiser et al.lll987t 
IGiovannini et al]l!993f) . We cannot directly compare our 
results with the projected angular distribution of radio 
emission resulting from the merger event, because our nu- 
merical code does not calculate the spatial dependence of 
the radiation emitted by particles. A merger model may 
be expected to produce softer radio emission at large 
distances from the cluster center, provided that the con- 
tribution of electrons and positrons formed in secondary 
nuclea r production from p rotons and ion interactions is 
small (|Miniati et al.ll200lfK This is because the speed of 
the merger shock varies according to the density contrast 
between the dominant and smaller merging cluster. Be- 
cause the core and maximum radii, r c and r max , of the 
minor cluster are less than that of the dominant cluster 
(Table 1), the speed of the shocked fluid becomes mono- 
tonically smaller at larger distances from the collision 
axis for a planar surface, essentially implying that a bow 
shock structure is created (see Ap pendix) . T his behav- 
ior is already apparent in Fig. 2 of ( BD 2003), where we 
found that shock strength is weaker at larger distances 
from cluster center, so that particles are injected with 
softer spectra in the outer parts of the merger, and even 
more so if they are away from the merger axis. Hydrody- 
namic simulations are needed to calculate the injection 
spectra at the same observing times, but it is clear from 
our results that softer emission will be injected at larger 
impact parameters while the merger event is in progress. 

In the case of Coma, excess gas located at a projected 
radial distance of ~ 0.5 Mpc from the cluster center, 
if interpreted as gas stripped from the merging cluster 
and left at the insertion point on the cluster periphery 
(R w 2 Mpc) , indicates that the angle between the plane 
of the sky and the orbital axis of the merger event is 
< 15 ° ijVikhlinin. Forman. k. Joneslll997t lArnaud et"afl 
120011) . Thus the propagation of the merger shocks are 
nearly along the line of sight, and the observed nonther- 
mal emission results from the line-of-sight integration of 
the particle emissivities. The lower density in the pe- 
riphery of the minor cluster as it merges with the larger 
cluster produces, in a planar geometry, stronger reverse 
shock emission with potentially harder injection at spec- 
tra large radii. If the two clusters were merging in a direc- 
tion transverse to the line-of-sight, then we would expect 
asymmetric emission related to the m erger geometry, as 
found in off-axis cluster-merger events l)Ricker fc Sarazinl 
2001). 

A bow shock structure, though with a rather low Mach 
numb er (M = 2.1), is seen in the merging galaxy cluster 
A520 i|Markevitch et al.ll2005l) that displays a prominent 
radio halo. Another is found in IE 065 7-56, which also 
shows a radio halo l|Govoni et~alll200l . This demon- 
strates that ongoing mergers are effective at accelerating 
primary electrons, though additional processes, such as 
turbulence, may contribute to acceleration in the radio 
halo. 
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In the picture studied here, the appearance of Coma is 
a result of an on-axis minor merger event where observa- 
tions are taken shortly after the passage of the centers of 
mass of the two merging clusters. This collision epoch is 
in fact indicated by the optical data, but could raise ques- 
tions of whether observations at this epoch are the result 
of fine-tuning. It seems likely that the merger shock is 
rapidly quenched after the passage and coalescence of the 
dark matter halos, so that radio emission is most likely 
seen precisely at the epoch at which we find the Coma 
cluster. This interpretation would also explain the rel- 
ative rarity of radio halos in clusters, which would be 
primarily limited to those which have undergone recent 
minor mergers. 

The symmetrical and unpolarized radio emissions of 
most cluster halos may be contrary to expectations from 
a merger event, which might be expected to be highly 
asymmetric and polarized as for radio relics. But for the 
minor mergers, which are the only ones with sufficient 
Mach speeds to make radio halos, a bow-shock "jellyfish" 
structure could be created that might be confused with 
a quasi-symmetrical radio halo. Viewing through both 
sides of the jellyfish could reduce the polarization. More- 
over, high resolution mapping of the structure of a radio 
halo could be unfolded to reveal the extended gaseous 
and possibly dark matter distributions. This prospect 
becomes more attractive with the development of high 
resolution, low frequency radio arrays to map the galaxy 
cluster radio halo emission. 

The difficulty for the cluster merger model is to avoid 
producing asymmetric radio emissions, especially for 
those events observed transverse to the collision axis. 
This problem cannot be alleviated by appealing to MHD 
turbulence that moves with the Alfven speed 

va = —j^= — 7 7 km s _1 (5) 



3 n_ 3 cm 3 



for 0.1-B_7/iG mean magnetic fields and 10 
mean proton densities, which is much slower than the 
thermal speed. Unless particles are transported without 
streaming losses, in which case they do not deposit en- 
ergy, energy is transported to the surrounding regions of 
the merger shock through MHD processes on timescales 
much longer than the merging timescale. A delayed re- 
acceleration model can, in principle, remedy the mor- 
phology problem. 

To m odel the radio ha l o em is sion from clus t ers o f 
galaxies. ISchlickeiser et al.l ill 9871) , IBrunetti et aTl 1120011) 
and iBrunetti et alJ l|2004) consider reacceleration of 
electrons via turbulence generated in a post-merger 
event. This approach may be supported observation- 
ally by detection of t urbulence in the core of Coma 
ijSchuecker et al.l 12004) . The total energy found in a 
population of relativistic nonthermal had rons is expected 
to be ~ 5% of the thermal energy pool l|Gabici fc Blasil 
120031 lBDll2003j) . inasmuch as most clusters are expected 
to experience several merger events in their formation 
history and protons with energies < 10 16 eV are con- 
fined within a c l uster of galaxies over a Hubble time 
ijVolk et alJll996UBerezinskv et alJll997|) . This nonther- 
mal pool of streaming protons could introduce turbulence 
that would accelerate electrons long after the merger 
events. Delayed reacceleration models are highly sensi- 
tive to the energy deposited into the nonthermal protons, 



however, so that if this energy content is unreasonably 
large, an episodic injection scenario, such as t he merger- 
shock model considered here, may be favored ijPetrosianl 
2001). 

The morphology question and the discrepant magnetic 
field values from Faraday rotation measurements com- 
pared to Comptonized CMBR/X-ray emission mode l val- 
ues of B = 0.2/iG l|Masil l200lt iPetrosianl l200l|) can 
in principle be reco nciled with the delayed turbulence 
reacceleration model l)Schlickeiser et aiTl 987: Blas l200H 
IBrunetti et al.ll2001|) . This model has its own difficul- 
ties, however. Second-order Fermi processes accelerate — 
and very slowly in the cluster merger environment — hard 
particle spectra. A pre-injection spectrum is required, 
such as pre-existing particles in the ICM. Particle spectra 
with the appropriate spectral index can explain the radio 
emission of the Com a cluster in a reacceleration model 
ijBrunetti et al J l2004j) . but have not yet been shown to 
explain the multiwavelength spectrum, including the ra- 
dio spectral softening with radius, by the same model. 
Pre-injection of nonther mal particles from AGN sources 
ijGiovannini et alj fl993) would also introduce asymme- 
tries that would have to be "washed out" in the tur- 
buelent reacceleration model. Indications that the ap- 
pearance of radio halos are cor related with recent cluster 
merger behavior l)Buotell2001l) are, however, apparently 
consistent with expectatio ns of the turbulence model 
ijCassano fe Brunet ti 2005) given assumptions about the 
energy density of turbulence with cluster mass. More 
detailed studies correlating peaks in radio halo emission 
and spectra with merging activity could distinguish be- 
tween a cluster merger and delayed reacceleration model 
but, most likely, both merger shocks and turbulent reac- 
celeration make important contributions to nonthermal 
activity in clusters of galaxies. 

Although the halo radio emission from Coma can be 
explained in the context of the cluster merger model, we 
do not claim that all radio halos can be accounted for by 
cluster merger shocks. A striking example that radio ha- 
los may have multiple origins is the recent XMM-Newton 
observation of the galaxy clusters A399 and A401. Both 
clusters have radio halos, but only the halo found in A399 
appears to be associate d with a shock front res ulting from 
a cluster merger event l|Sakelliou fe Ponmanl l2004). The 
origin of nonthermal electrons producing the radio halo 
associated with A401 could have a different acceleration 
mechanism, for ex ample, arising from accretion shocks 
( Enfil in et alJH99 8|) or revival of rel ic non-thermal parti- 
cles i|Enfflin fc Gopal-Krishnall200H) . 

Of central importance is to establish the strength and 
spectrum of the hard X-ray power-law component re- 
ported with BeppoS 'AX and RXTE, which normalizes the 
primary nonthermal electron power. Until these reports 
are confirmed, which may have to wait for the upcoming 
launch of Astro E-2, predictions for 7-ray detection of 
clusters of galaxies will remain highly uncertain. More- 
over, the relative efficiency to accelerate protons and elec- 
trons and the unknown fraction of residual nonthermal 
hadronic energy left from previous mergers add addi- 
tional uncertainties to predictions for 7-ray detection. 
First-order Fermi acceleration is, however, expected to 
accelerate protons with greater efficiency t han electrons 
becau se of the larger proton gyro-radii ijBaring et alJ 
1999), and past merger events can only add nonther- 



mal protons to the total. Thus our 7-ray predictions 
are expected to represent a lower range of the hadronic 
fluxes, provided that nonthermal hard X-ray emission 
has indeed been detected with BeppoSAX and RXTE. If 
correct, then we predict that the Coma cluster will be 
detected with w 5<r significance with GLAST (see Ap- 
pendix), and marginally or significantly detected with 
VERITAS and HESS. 

Our results can be comp ared with a recent theoretical 
study ijReimer et alJl2004|) which concludes that Coma 
will be detected with GLAST but not with the IACTs 
. The level of the high-energy radiation fluxes, which is 
made primarily by ir° decay emission, depends crucially 
on the spectr um and high-energy c utoff of the nonther- 
mal hadrons ijAtovan fe Volkl l2000j) . which we argue is 
correctly calculated here. Observations of 7-ray emission 
from Coma will reveal the relative efficiency of nonther- 
mal proton and electron acceleration. 

High resolution radio maps from low-frequency (< 100 
MHz) arrays will be important to determine the nature 
of the observed radio emission in Coma. In the clus- 
ter merger scenario, energetic protons will diffuse far- 
ther from their point of origin than lower energy protons, 
and these have had a much longer time than the merger 
timescale to diffuse outwards. At sufficiently large radii, 
there will be a transition from radio emission of primary 
to secondary origin. Thus the cluster merger model indi- 
cates that the spectral index of Coma's radio halo softens 
with radius until the radio emission is dominated by lep- 



ton secondaries from hadrons, where it would begin to 
harden. A hardening in the volume-integrated mm-radio 
spectrum of Coma is likewise expected from a hadronic 
component. 

The > TeV 7-rays from Coma will also be atten- 
uated by the diffuse extragalactic infrared radiation 
ijGabici fc B lasi 2004|). Our numerical results from clus- 
ter merger model indicate that acceleration of hadronic 
cosmic rays by cluster merger shocks will produce a 
very high-energy 7-ray emission component extending to 
> 10 16 eV. Compton cascading and reprocessing of the 
cluster halo ultra high-energy (VHE) 7-rays will, for cer- 
tain values of the intracluster magnetic fields, produce 
extended cascade 7-ray halos. Detection and mapping 
of such halos would be an important probe to measure 
ICM magnetic fields. 
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APPENDIX 

ANALYTIC CLUSTER MERGER PHYSICS 

Radiation Processes 

Suppose that the ICM of Coma is characterized as a quasi-spherical hot cloud of gas and plasma entraining magnetic 
field and nonthermal particles that were injected by shocks formed by merger events during structure formation. 
Assume an isotropic nonthermal electron pitch-angle distribution. On the size scale of the telescope's resolution, we 
also assume that the magnetic field direction is randomly ordered, though with a mean magnetic field intensity B that 
is roughly constant through the cluster. 

An electron with Lorentz factor 7 emits photons with frequency v syn (Hz) = 2.8 x 10 6 _B(G)7 2 . For magnetic fields 
measured in units of 10 -7 -B_7 G, and for radio measurements at z^ghz GHz, electrons with Lorentz factor 

7 S yn = 6xl0 4 ^|^, (Al) 

Electrons with Lorentz factor 

7c = 2.5 x io 4 yJj^— , (A2) 

Thomson-scatter CMBR photons to X-ray energies (e = 0.1 corresponds to a photon of « 50 keV). 

Let Pqhz and e denote the radio frequency and dimensionless hard X-ray photon energy where the spectral slopes 
of the nonthermal radiation are the same. For Coma, we find that a non thermal photon flux with p hoton number 
spectral index = —1.8 (a = 0.8) provides a good fit to the BeppoSAX l|Fusco-Femiano et al.l 120041) and (RXTE) 
(Rep haeli fc Gruberll2002l) observations (Fig. EJ. The radio index a = 0.8 at 30 MHz energies (Fig. 0}. The interpre- 
tation that the nonthermal radio emission is synchrotron radiation and that Coma displays an associated nonthermal 
hard X-ray component from Compton-scattered CMBR implies magnetic fields from the condition that the primary 
Compton and primary synchrotron fluxes should have the same spectral indices for electrons with the same energies, 
that is, 7s yn = 7c- Hence 

5_ 7 ^5.8 — (i + z ) . (A3) 
e 

For Coma, therefore, B^j = 2.3, in good agreement with the value B^j = 2.2 used in the numerical model. The 
discrepancy between this value and val ues derived from Faraday rotation measurments could be explained by mul- 
tiple zones (Blasi ll200lt|Pet rosianll2"001|) . or the attribution of the hard X-ray excess to nonthermal bremsstrahlung 
(jSarazin fc Kemrmeril200'ci() . Confirmation of the Beppo-SAX and RXTE results is especially crucial here, because the 



magnetic field can be larger and the electron intensity correspondingly smaller if there is no nonthermal X-ray tail to 
explain. 

The timescale for Thomson-scattered CMBR and synchrotron cooling can be written as 



2.5 x 10 6 



Myr. 



7 [(l + z) 4 + 10- 3 B2 7 ] 

The lifetime for electrons that emit synchrotron radio emission with frequency z^ghz GHz is 



tc s (7s; 



4VB-?/ ^GHz 



Myr. 



(1 + zf + 10- 3 B 2 _ 7 

The lifetime for electrons that emit hard X-ray photons with dimensionless energy e = 0.1e_i is 

320^(1 + z)A-i ^ 
^ (7c)= (l + z)4 + 10-3^ 7 Myr - 

These timescales are generally shorter than dynamical time scales, as we now show. 



(A4) 



(A5) 



(A6) 



Scaling Relations for Cluster Dynamics 

From elementary considerations, the gain in kinetic energy when the minor cluster, treated as a test particle in the 
mass distribution of the dominant cluster, falls from radius n to radius r 2 (< ri), is M<2,v\j2 = GMiMiiy^ — 
Thus v = yj 2GM\ I r2 when r2 < n, so that 

v 2 « 6000 



(r 2 /0.257 Mpc; 



km s 



(A7) 



where r 2 is equated with the core radius of the dominant cluster (see Table 1). The characteristic merger time t is 
determined by the acceleration a = GM\jr\ at the outer radius. Because T\ rj a 2 i 2 /2, 



t 



2r\ 



3/2 

660 ^E£ Myr, 



y GM\ M 15 

where Mi — 10 15 Mq and r\ = rM pc Mpc. The available energy in the collision is 

M 2 



c GM\M<i 8 x 10 63 

t f» f» M15 

r r Mpc 



10 14 M, 



ergs 



(A8) 



(A9) 



The long timescale of merger events compared to the timescale for electrons that emit GHz radio emission, eq. (| A5|> . 
means that these electrons had to have been accelerated recently compared to the duration of the cluster merger event. 
The sound speed 



lk B T 



1200 



k^T x 
lOkeV 



km s 



(A10) 



where the mean mass per particle is taken to be (m) = 0.6m p . From eq. l)A7j) . this implies shocks with moderate Mach 
numbers M. ~ 5. 

Analytical Cluster Dynamics 

We treat cluster dynamics analytically (see also lRicker fc Sarazinll200"ll) . as a check on the validity of the numerical 
treatment. The merging activity in the Coma cluster of galaxies is again treated in the test-particle approximation for 
the minor cluster with mass Mi, though now using the /3-model approximation to the isothermal mass distribution 
ijCavaliere fc Fusco-Femianolli976D given by 

21 "3/3/2 



P(r) = A) 

for the dominant cluster, with po = CPo,ii from eq. The mass of cluster 1 interior to radius r is 



, r 
1+ - 



Mi(< r) = 47r p r s c 



•/r c 



dx 



1 + X 2 



(All) 



(A12) 



We approximate (3 = 2/3 for Coma, which is close to its value, (3 = 0.705, determined from observations. The outer 



radius 



is normalized to the total mass M\ of the dominant cluster, giving 



Mi = Amrt 



arctan 



Po 



(A13) 
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Solving for ( — 21, corresponding to a ratio of normal mass to CDM mass of 5%, gives r max = 0.82 Mpc. When 
( = 11, for a 10% ratio of normal mass to CDM mass, r max = 1.3 Mpc. 

The change in potential energy of the test-particle cluster 2 when falling through the gravitational potential of cluster 
1 is 



1 



AU=-M 2 {Av) 



-4irGp M 2 r 3 c 



dr r 



r/r 



dy 



V 



( 1 + y 2)3/3/2 



provided r max >r\ > r 2 . For the case (3 — 2/3, eq. (|A14J) can be solved to give 



(Av) 2 = AwGpor 2 



In 



i + vi 



1 



vl 



( arctanj/! arctany 2 



(A14) 



(A15) 



V Vl V2 
where y 1(2) = r 1(2) /r c . 

Taking y 2 — * 0, we obtain Aw = 4760 and Av = 3750 km s _1 for £ = 21 and 11, respectively. A sound speed c s = 1100 
km s _1 of the dominant cluster (Table 1 and eq. jAlOp implies Mach numbers M. = 4.37 and 3.44, compression ratios 
X = 3.5 and 3.2, and injection indices s = 2.2 and 2.4 for ( — 21 and 11, respectively, using the standard formulas for 
nonrelativistic shock acceleration in the test particle limit. The injection index s = 2.2 is in good agreement with the 
numerical results (see Fig.^l. 

We instead consider the dark-matter density profile distribution given by a 

-0 



p{r) 



1 



(A16) 



where r s is the scaling radius and p s is the characteristic overdensity. For the NFW ijNavarro. Frenk. fc Whitelll997j) 
model, (a, /3, 7) = (1, 2, 1), and the mass of cluster 1 interior to radius r becomes 



Mi(< r) = 4:irp s r 3 s \ dx 
Afx(<r ; 







2 = Airp s r^ In 
) at the outer radius r : 



1 + — I - 



r/r s 



1 + r/r, 
inn,, one obtains that r 



(A17) 
843 kpc for Mi 



By normalizing to the total mass Mi 

from Table 1, with p s = 1.13 x 10~ 25 g cm~ 3 (using r s — 459 kpc and c = 5.42 given bv lEttori. De Grandi. fc Molend) 
2002). Deriving the speed following the approach leading to eq. I|A15|I gives 



(Av) 2 = 8TTGp s r 2 



ln(l + y) 



y 



ri/r B 



T<2 1 r s 



(A18) 



In the limit r 2 <C r s and r\ = r max , one finds that Av = 4040 km s , implying a Mach number of M. = 3.7, 
intermediate to the values found for the /3-model profile. Thus, our conclusions should not be much different whether 
we consider a /3-model or NFW profile. 



In B 



Shock Structure in Merging Clusters of Galaxies 
l)2003|) . we derived the Mach numbers Mi and M2 of the forward and reverse shocks, respectively, given by 

2 v 



and 



3 c 2 




(A19) 



(A20) 



Here ci and c 2 are the sound speeds in the dominant and merging clusters, respectively. The value of v is calculated 
by iteratively solving 

Ml ni(t) _ 1 + 3M± 2 (t) fv (t) - 



, , . (A21) 

P2n 2 (t) l + 3Mz 2 (t) 1 

using eqs. (|A19|) and l|A20|) . The terms ni(t) and n 2 (t) are the ICM number densities given by eq. IjAllfl . and pi and 
P2 are the mean masses per particle in the dominant and merging clusters, respectively. 

A core radius of 0.15 Mpc is used for the minor cluster in the calculations. Becasue the core radius of the minor 
cluster is smaller than that of the dominant cluster, the relative densities along the shock front decrease away from 
the collision axis joining the two merging clusters. The shocked fluid speed also monotonically decreases away from 
the collision axis, as can be seen by numerically solving the above equations. This produces a weaker forward shock 
and a bow shock structure that will inject softer electron spectra at greater angular distances from the center of the 
Coma cluster. 

The speed v of the shock ed fl uid and the Mach numbers of the forward and reverse shocks calculated from these 
equations are shown in Fig. IA1I as a function of density contrast ni/n 2 in a planar geometry. The relative speed Vq 
of the two clusters is set equal to 6000 km s -1 . The gas temperatures of the two clusters are given in Table 1. When 
ni/n 2 is larger, the speed of the shocked fluid is smaller. Consequently a bow-shock structure is formed. 
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Fig. Al. — The ratio v/vo of the shocked fluid speed to the merger speed i>o, and the Mach numbers A4i and M2 of the forward and 
reverse shocks, respectively, are plotted as a function of the ratio of the densities of the dominant and merging clusters. 

Detection of Coma with GLAST 

Using the expressions derived in Section 4 of IBDI l)2003j) . we estimate the detection significance and number of 
photons that GLAST will detect from the Coma cluster during t yr years of GLAST observations in the scanning 
mode, using an exposure factor for Coma of 20%. 

Figs. |21 and show that the 7-ray flux spectrum from Coma can be expressed as 

" i - = 1 ° 43t «(ioow)'""^ s "- < A22 » 

where a v — —0.2 and £43 < 0.3 (in the best case still allowed by the high-frequency radio data, r/ p = 20% and 
L 43 «0.6). 

The number of photons with energies > E 1 detected from the Coma cluster with GLAST is 

N s (> E y ) - 35 Vl^GeV)]- 1 - 04 . (A23 ) 

More photons are detected at lower energies, but this expression is only valid above the ir° —> 2j bump, which is found 
at a few hundred MeV in a vF v representation. 
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The best detection significance is found by assuming that the emission from Coma is point-like, which holds until 
the point-spread functi on of GLAST is smalle r than the angular extent of Coma. The Coma C radio source has an 
angular extent of w 20' ijGiovannini et al1l!993f) . which is comparable to the psf of ss 2 GeV photons. The n a detection 
significance for Coma is found to be 

n a = 5.4 (ff ) [^(GeV)] 01 . (A24) 

These two expressions show that Coma will be significantly detected with GLAST if the energy injected in nonthermal 
hadrons is ~ 10 x the energy in nonthermal electrons (or if there is significant energy in nonthermal protons left over 
from previous merger events). The number of detected photons, which may be as large as a few hundred, would permit 
spectral analysis to determine if the spectrum is consistent with a —2.2 photon spectral index at GeV energies. If the 
nonthermal hadron energy content has been overestimated, however, it remains possible that that GLAST will only 
weakly detect the Coma cluster of galaxies. 
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